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ABSTRACT 

We present a panoramic atlas of SpitzerfiRSi spectra of extragalactic sources collected 
from the recent literature, with value added measurements of their spectral features obtained 
in a homogeneous and concise manner. The atlas covers the full spectrum of the extragalactic 
universe and includes star forming galaxies, obscured and unobscured active galaxies, lu- 
minous and ultra-luminous infrared galaxies, and hybrid objects. Measured features such as 
the polycyclic aromatic hydrocarbons, the strength of the silicates in emission or absorption 
around 9.7 [xm, rest-frame monochromatic luminosities or colours, combined with measure- 
ments derived from spectral decomposition, are used to establish diagnostics that allow for 
classification of sources, based on their infrared properties alone. Average templates of the 
various classes are also derived. The full atlas with the value added measurements and ancil- 
lary archival data are publicly available at http ://w w w.denebola .org/ atlas , with fuU references 
to the original data. 
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1 INTRODUCTION 

The mid-to-far infrared spectral range holds the key to the physi- 
cal processes that shape the life and properties of galaxies, namely 
star formation and nuclear activity. This spectral domain, and with 
the exception of the very bright emission lines of the brightest 
sources, was almost unaccessible from the ground and until 15 
years ago, when the sensitivity and large wavelength coverage of 
the Infrared Space Observatory (ISO) allowed the detailed study 
of large numbers of nearby starburst galaxies and active galactic 
nuclei (AGN) and lead to the first mid-infrared (MIR) diagnos- 
tic tools. ISO studies observed, among other things, the presence 
of polycyclic aromatic h ydrocarbon (PAH) features in the majority 



or polycyclic aromatic n ydrocarbon (tyvhl ) features in tne ma|ority 
of star-forming galaxies dLutz et alJll998l : rRigopoulou et al.l I999I) 
and their more scarce appearance in AGN (e.g. lsturm et alj|2000h . 
as wel l as the presence o f a silicate feature in absorption in type 2 
AGN dClavel et al.ll20 00'). 

Though the studies carried out with ISO brought a break 
through in our understanding of star formation in galaxies, it was 
the Spitzer Space Telescope with its MIR cameras IRAC and 
MIPS, and the InfraRed Spectrograph (IRS) that revolutionised our 
knowledge on the starburst and AGN phenomena and their con- 
nection. IRAC broad band MIR photometry alone allowed for a 
rough classification among AGN and a distinction between active 



e-mail: ahc@deneboIa.org 



and star-forming galaxies (e.g. Lacy et al]|2004l : IStern et al]|2005l : 
iHatziminaoglou et alj20"o^l2009l) . while IRS observations of hun- 
dreds of extragal actic sources provid ed with more elaborate diag- 
nostic tools (e.g. ISpoonetaljboOTi) and showed that their MIR 
spectra are richer and more complex than once assumed. The new 
discoveries, however, raised new issues that models are now chal- 
lenged to reproduce. The silicate emission feature in type 1 AGN 
peaks at longer wavelengths than the 9.7 |xm predict ed by the mod- 
els (e.g. lSiebenmorgen et aI.ll2005l : lHao et al.ll2()()^: silicate emis- 
sion has also been observed in type 2 AGN (Sturm et al. 200^; 
iTepIitz et aLll2006l : iNikutta et al.i ,2009) in apparent contradiction 
to at least some of the suggested dust geometries; and PAH e mis- 
sion is present in many quasars (e.g. ISchweitzer et al.ir200(j ). de- 
spite previous theories postulating their destruction in the vicinity 
of the active nucleus. Meanwhile, accurate measurements of PAH 
features in star-forming galaxies are made difficult by the presence 
of other features, such as the deep silicate absorption at 9.7 |,im or 
an underlying AGN continuum. 

Understanding the processes linking different galaxy popula- 
tions and finding their distinctive characteristics requires the com- 
parison of the properties of large samples, selected with different 
but well understood criteria. Studies of individual samples usu- 
ally pursue specific goals, and basic MIR properties for diagnostics 
(such as the strength of the silicate features and PAH bands) are 
measured in ways adapted to the needs of each individual study, but 



2 Herndn-Caballero 6^ Hatziminaoglou 



not necessarily consistently across samples, making a direct com- 
parison of results painful, if not impossible. 

With this in mind, we compiled a list of 739 SpitzerfiRS spec- 
tra of extragalactic sources, the largest of its kind to date, previously 
scattered in the recent literature. The sample, though not statistical, 
comprises a variety of objects, namely type 1, 2 and intermediate 
AGN, ULIRGs, starburst and sub-mm galaxies, with a large span 
of spectral properties such as PAH and silicate features, continua 
and distinct physical characteristics. The MIR properties were me- 
assured in a concise way for all objects and can thus be directly 
compared. The paper is structured as follows. §2 describes the sam- 
ple selection, overall properties and collection of ancillary data; §3 
presents the extraction of the IRS spectra and calibration proce- 
dure followed; sections 4 and 5 detail the methods used to measure 
MIR spectral properties and the composite spectra obtained for the 
various subsamples, respectively; §6 reproduces a few common di- 
agnostic diagrams and finally §7 summarizes our conclusions. 

The full collection of IRS spectra with references to the orig- 
inal data, ancillary data and additional measurements are available 
at http://www.denebola.org/atlas 
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Figure 1. Redshift histogram for the 695 sources with spectroscopic red- 
shifts in the sample. The red bars indicate optical redshifts, while the green 
bars represent the MIR ones. Bars for the first two redshift bins are truncated 
to improve readability. 



2 SAMPLE SELECTION 

The sample hereby presented is drawn from papers in the 
literature presenting SpitzerfiRS spectroscopy of active and 
star-forming galaxies, and contains a wide variety of tar- 
gets including Starbursts, Seyferts, LINERs, ULIRGs, QSOs, 
SMGs and radio galaxies, a mong others . A tot al 739 unique 
source s are selected f r om Brand et al. I ( 12008 ). ISrandl et al 
20061) , iBuchanan et alj (2006h. ICao et al .l Hooi), lOasvra et al 



2009l),lDeoetalJ l l2007h 



iHaoetal.' ('2005), 'Heman-Caballero et al. (2009), 'Hiner et al 
( ^2009 ), Imanishi et al. (2007, 2010), Imanishi (2009), Lacy et al, 
20071), iLeipski et al.l (l2009h. iLutz et al. (2008^, i Maiolino et al 



2005), 



Farrah et al 



:tal.l (I20U«I ). lUasyra et al 
]2009), Haa s et alj (1200^ 



rr 



Shi et al.' ('2006*), 
(2005, 2006), Tommasin et al 



20071), iMenendez-Delmestre et al.h2009 ). Murph v et al.l i 200S ) 
Netzer et al.l j2007l) . IPolletta et alj jjo OS), Pope et al.1 120081) 

Siebenmorgen et al. (2005), Sturm etal 
(2008), VaUanteet^ 



Weedman et alj 1 200: 



Wu et al4 i2009i . lYan et al.l ^20071) and lZakamska et al.l j2008h" 



2006 jfbl 



IWillett et al. 



Table[T]summarizes the type and number of spectra and obser- 
vation mode for each subsample. For a more detailed description of 
the selection criteria and properties of the individual subsamples we 
defer the reader to the corresponding papers. 

SpitzerjlRS data for sources in the sample include observa- 
tions with both the Low-Resolution (R~100) and High-Resolution 
(R~600) IRS modu les, all performed in s taring -mode, ex c ept fo r 
the sources from iBuchanan et alj j2006l) and IWuetalJ j2009l) . 
where observations were performed in spectral mapping mode, but 
only the spectrum obtained with the slit placed in the galactic nu- 
cleus is given. Maximum spectral coverage is 5-38|xm in the Low- 
Res module and 10-35|^im in the High-Res one, but in many faint 
targets only a subset of Low-Res modules is applied, resulting in a 
reduced spectral coverage that is not uniform even within the same 
subsample. 

Sources at intermediate and high redshift appear unresolved 
to the IRS, and thus the spectra integrate the emission of the whole 
galaxy. By contrast, in nearby targets the spectra only sample the 
nuclear region. Our selection does not put constraints on redshift, 
not even the existence of a consistent redshift estimate, but individ- 



ual subsamples usually target certain redshift range, or indirectly 
bias the selection through flux or colour constraints. 

Spectroscopic redshifts (either from optical or MIR spec- 
troscopy) are available for 695 of the 739 sources. Redshifts are 
obtained from the literature and cross-checked with the NASA Ex- 
tragalactic Database (NED). 525 redshifts come from optical spec- 
tra, while 170 sources (mostly obscured high-z targets) rely on a 
redshift estimate from their IRS spectrum. In targets where both 
optical and MIR spectroscopic redshifts are available, priority is 
given to the former, except in cases of strong discrepancy in which 
the optical redshift is flagged as unreliable or the MIR spectrum 
shows very clear PAH bands. The remaining 44 sources have no 
published spectroscopic redshift, or it is uncertain. They are not 
removed from the sample, but will be ignored in the subsequent 
analysis. 

About half of the sample, 347 sources, have redshifts below 
z=0.3, while the rest spans a large redshift range up to z~3.7 (Fig- 
ure [TJ. Spectral coverage in rest frame wavelength is in all cases 
within 1.2 and 38 |,im, with 389 sources sampling the 5.5-14 |xm 
range in its entirety and 550 the critical 8-I2|.im range, which con- 
tains both the 10 [im silicate feature and 11.3 |xm PAH band (see 
Figure |2](. 

An optical spectroscopic classification is also available for 
most sources. A positional match with the 1 3 th version of the Veron 
Catalogue of Quasars &. AGN (VERONCAT; IVeron & Verori2006l) 
yields 361 matches; for these sources, we adopt the activity type 
designation in the VERONCAT. For the remaining ones we adopt 
the optical classification in NED, when available, or at the very 
least a broad class (e.g.: ULIRG, SMG) either from NED or the 
individual papers the spectra are taken from. Some sources are put 
in more than one classes. For example, a starburst galaxy with a 
Seyfert 2 nucleus will retain both starburst and Seyfert 2 clasifica- 
tions. 

To deal with the large variety of spectral types, we group them 
into 10 categories, namely: Seyfert 1 (Syl), Seyfert 2 (Sy2), in- 
termediate Seyfert type (Syl.X), LINER, Quasar (QSO), type 2 
Quasar (QS02), Starburst, ULIRG, SMG and "other" when none 
of the above classifications apply. 
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Figure 2. Restframe wavelength coverage for the 695 sources with spectro- 
scopic redshift in the sample. Each point represents one source, with the X 
and Y coordinates indicating the shortest and longest wavelengths sampled 
in the IRS spectrum, respectively. Symbols indicate spectral coverage in 
the 5-15(i.m restframe range: >95% (blue circles), between 60% and 95% 
(green boxes) and <60% (red triangles). 



ertheless, the uncertainty that this introduces to the wavelength cal- 
ibration is about an order of magnitude smaller than the spectral 
resolution (R~100) in the Low-Resolution module of IRS, and thus 
its impact is negligible. 

Since the uncertainties in the flux density of the spectra are not 
shown in most of the figures, we devised a procedure to estimate 
them based on individual exposure times, obtained from queries 
to the Spitzer archive when they were not indicated in the corre- 
sponding papers, using the SPEC-PET exposure time calculator. 
Flux uncertainties estimated by SPEC -PET are corrected for wave- 
length sampling (SPEC -PET assumes R=50), with an additional 
V2 factor to account for the noise introduced by standard back- 
ground subtraction. We use two different methods to evaluate the 
accuracy of these flux uncertainty estimates: i) compare with the 
known uncertain ties in the spectra from papers where this informa - 
tion is available: 'Hernan-Caballero et al (iQQ't). ' Yan et alj (l2007h , 
[iVIenendez-Delmestre et al. (2009) and Murphy et al.l ( |2009|) : li) 
compare with the RMS noise of the c ontinuum i n spectra that are 
plotted without smoothing (e.g. [Weedman et al] 1200^ : iLutz et al] 
I2OO8I : IPope et alj|2008l ; IZakamska et al.ll2008h . 

Both approaches indicate that our flux uncertainty estimates 
are correct within a factor 2 in most sources. Even if a higher ac- 
curacy would be desirable, the statistical results drawn from the 
sample are little, if at all, afl'ected. For an in-depth analysis of in- 
dividual targets authors should use instead the latest version of the 
data products available at the Spitzer archive. 



We set a threshold value of vL^(7[,im) = 10"** erg s"' (roughly 
the output of the least luminous ULIRG of this sample), above 
which all type 1 and type 2 AGN are considered as QSO or QS02, 
while keeping the Seyfert 1 and Seyfert 2 designations for the lower 
luminosity Seyferts. The LINER class does not include LINER 
cores of ultraluminous and starburst galaxies, as their MIR spec- 
trum is typically dominated by the host galaxy. 



4 MEASUREMENTS 

We provide measurements for several MIR spectral features that 
are commonly used to estimate the star-formation or nuclear activ- 
ity, their relative contribution to the energy output of extragalactic 
sources or the amount of obscuration; namely PAH features, the 10 
^im silicate feature and rest-frame monochromatic luminosities. 



3 EXTRACTION OF THE SPECTRA 

To obtain the IRS spectra we have not gone through the normal 
process of querying the Spitzer archive, reduce the raw Spitzer/IRS 
data and extract the spectra. Instead, we have relied on Postscript 
figures contained in papers submitted to the arxiv.org preprint 
archive. This approach has several shortcomings that we describe 
later, and is not optimal for a detailed analysis of individual tar- 
gets. Nevertheless, it provides a cost-efficient means for performing 
statistical studies on the largest starburst and active galaxy sample 
compiled to date. 

Sets of points representing the spectra in the figures are trans- 
formed from Postscript coordinates to wavelength and flux values 
using software developed by the authors. Paper figures represent 
the spectra in many difl'erent fashions: the independent variable is 
usually wavelength, whether in the observed or rest frame, but it 
can also be frequency. For the dependent variable there are several 
values of choice: f^, f^, yf,,, vLy, etc., that in addition are expressed 
in different units systems and plotted in linear or logarithmic scale. 

To homogenize the sample, we convert all spectra back to 
the same units provided by the Spitzer IRS reduction pipeline: ob- 
served wavelength in micron and spectral flux density in Jy. The 
accuracy with which the original wavelength and flux values are 
recovered is limited by the resolution of the Postscript figure; nev- 



4.1 PAH features 

PAH emission is a key tracer of star formation, with the lu- 
minosity of the PAH bands frequently used as an indirect es- 
timat o r of the Star Form a tion Rate (SFR; e. g. Rigopoulou et al] 
1 19991 ; |PeetersetalJl2004l ; iBrandl et alj |2006| ; iFarrah et al] I2OO8I) 
while their ratios relate to the ionization field in the interstel- 
lar medium jPraine & 0120011 ; iRapacioli etal]|2005l ; iBrandl et all 
I2OO6I ; iGalliano et ai]|2008l) . For nearby galaxies with strong PAH 
bands and high signal-to-noise (S/N) spectra, it is possible to per- 
form a detailed modelling of th e continuum and Lorentzian p rofiles 
of the PAH bands (e.g. ISmith et al]|2007l ; iTrever et alJlIol O). For 
higher redshift sources, a similar analysis can be performed, pro- 
vided adaptatio ns are made to account for AGN emission and dus t 
extinction (e.g. ISaiina et al]|2007l ; iHernan-Caballero et alj|2009l) . 
but as the equivalent width (EW) of the feature, the S/N of the spec- 
trum and the wavelength coverage decrease, degeneracy among so- 
lutions quickly becomes a major source of uncertainty. Simpler 
methods such as spline or linear interpolatio n of the local con- 
tinuum and integration in the PAH band (e.g. Genzel et al.' 199^; 
Rigopoulou et al. 1999; Brandl et al. 2006; Spoon et al. 2007) are 
much more resistant to weak or noisy features, but significantly 
underestimate PAH luminosities due to losses at the wings of 
the Lorentzian profiles and overlap b etween adjacent PAH bands 
JSmith et al.ll2007l; iTrever et alj|2010l) . 
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As already mentioned, the basic idea beiiind tiiis work is to ap- 
ply the same procedure when measuring the various features across 
all samples. We therefore opted to use an intermediate approach, in 
which the continuum is locally interpolated and the feature flux is 
integrated in a narrow band, but both measurements are later cor- 
rected for the expected width and shape of the feature's profile. 

The measurement proceeds as follows: we select two narrow, 
continuum bands at both sides of each PAH feature, and estimate 
the average flux, f^, in these bands. We perform a linear interpola- 
tion to estimate the continuum underlying the feature, and subtract 
it from the spectrum. Finally we integrate the residual in a band 
centered at the expected wavelength of the peak of the PAH feature 
to obtain the integrated PAH flux. In order to maximize the S/N in 
the integrated PAH flux and to reduce the uncertainty in the under- 
lying continuum, we use a rather narrow integration band (which 
loses a significant fraction of the PAH flux in the wings of the pro- 
file) and place the continuum bands in its close vicinity, without 
trying to avoid contamination by the wings of the PAH feature (see 
parameter values in Table [3}. 

The bias introduced by this procedure is compensated by 
a correction factor on the assumption that the PAH profile is 
Lorentzian with a known full width half maximum (FWHM). We 
estimate the error that this procedure introduces by meassuring 
line and continuum fluxes in a set of model continuum+Lorentzian 
spectra covering a wide range of EW, FWHM, and S/N values. We 
find that a 10% increase in the FWHM causes a 5-6% drop in the 
PAH flux, with no observable dependency on the EW or S/N of the 
spectrum. 

Uncertainties in the PAH flux and the underlying continuum 
for each source are estimated by performing Monte Carlo simula- 
tions. Finally, we estimate EWs of the PAH features dividing the 
integrated PAH flux by the interpolated continuum at the center of 
the feature, and convert PAH fluxes to luminosities assuming a con- 
cordance ACDM cosmology (A = 0.7, Q. = 0.3 and Ho = 70 km 
s"' Mpc"'). 

Since this approach is not valid for the overlapping 7.7 and 
8.6 PAH bands, we perform measurements only for the 6.2 and 
11.3 i-im PAH features (which incidentally are the most used for 
diagnostics). The 6.2 |,im feature is measured in 587 sources, that 
at 11.3 [im in 577, and both of them together in 481 sources. 

4.2 Silicates 

The 10 |,im silicate feature, which can appear in emission or 
absorption, provides insight into the geometry of the dust dis- 
tribution along the line of sight of the star-forming regions 
or t he active nucleus as wel l as the amount of obscuration 
(e.g. Ilmanishi & Malonevll2003l ; Isturm et al. 2005; Imanishi etjal] 
12007 *). As with the PAH bands, the main difficulty in estimating 
the strength of the 10 |xm silicate feature is the identification of the 
underlying continuum. In sources with weak or no PAH emission, 
we interpolate linearly in log(/l), log(f^) between anchor points at 
restframe 8.1 and 14.0 |,im; but whenever there is significant PAH 
emission (rpQR > 0. 1, see we substitute the first anchor point 
by 5.55 |xm in order to avoid the distortion introduced by the wings 
of the 7.7 and 8.6 |xm PAH features. 

We subtract from the spectrum the interpolated continuum, 
and integrate the residual (silicate feature) in the 9-11 |,im range. 
The sign of this quantity allows to determine whether the fea- 
ture is in emission (integral > 0) or in absorption (integral < 0). 
The theoretical wavelength of the silicate peak is 9.7 |xm, but it is 
known that in many quasars and Seyfert galaxies the feature peak 
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Figure 3. Comparison between two estimates of the silicate strength: 5sii, 
obtained by the standard procedure (see text) and C7-10-14, calculated as a 
three-band colour index. Sources with low uncertainties in the estimation of 
5sii (eSsii < 0.2) are shown in filled circles, while open circles represent the 
remaining objects. 



is shifted to longer w avelengths when it appe ars in emission (e.g. 
ISiebenmorgen et aP 12005 : Netzer et aDl2007h . In addition, many 
star-forming sources show strong emission in the S(3) 9.665 [xm 
transition of molecular hydrogen, which contaminates the peak of 
the silicate absorption profile. 

Since many spectra are too noisy to accurately measure where 
the peak of the silicate feature actually occurs, we simply assume 
it is located at A,, = 9.8 [xm for sources with silicate absorption and 
Ap = 10.5 ^m for silicate emission spectra. A detailed analysis of 
the 10 and 18 ^m silicate features in the AGN-dominated sources of 
this sample, including measurement of the peak's wavelength, is the 
topic of a specialised study (Heman-Caballero et al. in preparation). 

The feature's peak flux density is used to calculate the silicate 
strength, 5sii = \n[F(Ap)fFc{Ap)]. As in the PAH bands, uncertainties 
in the silicate parameters are estimated using Monte-Carlo simula- 
tions. 5sii measurements are obtained for 537 sources. 

We also define the "three-band" colour index, Cv-io-m, as: 

C7-10-14 = -2.5 tog(-^^^) 

\ Jio I 

This parameter is an alternative, easily obtained estimator of 
the silicate strength, that shows a good correlation with 5sii in 
sources of all types (Figure[3]l. 

In sources with little or no PAH emission, 5sii and Cv-io-m are 
mostly interchangeable, while sources with strong PAH emission 
show a small but consistent bias towards Cy-io-u < S^^. This may 
be produced by a steeper slope of the starburst continuum compo- 
nent shortwards of ~8.5 |xm that decreases the interpolated 9.8 |,im 
continuum for S^w (and increases thus the S^w value), and also from 
contamination by the wings of the 6.2 and 7.7 [j,m PAH features 
in the 7.0 [rm flux, that reduces the Cv-io-m estimate. The differ- 
ence, Asii = ^sii -Cv-io-w increases with the strength of the starburst 
component as measured by EW62> and can be as high as ~ 0.5 in 
extreme cases. 
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4.3 Monochromatic luminosities 

Rest-frame monochromatic luminosities, vLv, allow for the com- 
parison of spectral energy distributions of sources at diiTerent red- 
shift without the need for K-corrections. We estimate them based 
on the IRS spectra at rest-frame wavelengths: 2.2, 3.6, 5.6, 7.0, 9.0, 
10.0, 12.0, 14.0, 15.0, 18.0, 20.0 and 25.0 ^m. To this end, each rest 
frame spectrum is integrated in a narrow band (between 0.5 and 2.0 
|,im wide depending on the wavelength) centered at the nominal 
wavelength, and divided by the band width. Table |2] summarizes 
these measurements. 

Rest frame MIR colours and colour indices are also calculated 
and tabulated. Such quantities give a quick but accurate idea of 
the shape of a spectrum. The [5.6-15] |,im colour, for instance, is 
a proxy for MIR continuum slope, while the three-band 7-10-14 
|,im colour index (Cv-io-m) can be used as a substitute estimate of 
the silicate strength in noisy or poorly sampled spectra (see 34. 2t . 



5 SPECTRAL DECOMPOSITION 

Decomposition of MIR spectra in several spectral components is a 
powerful tool that can provide considerable insight into t he physics 
of the sources (e.g. lLaurent et alj2000l : lTran et alj200lh . 

MIR emission from active and star-forming galaxies arises 
mostly from HII regions, Photodissociation Regions (PDRs) and 
hot dust heated by energetic X-ray-to-UV photons from an AGN, 
so we have used a simple model comprising the superposition of 
three spectral templates (HII, FDR and AGN) obscured by a screen 
of dust. 

The spectral decomposition is performed by fitting the 5-15 
|xm restframe range of the spectrum to a parametrized F,i(A) of the 
form: 

FJA) = e-'"''''>[aJj,oNW + a2fmiW + «3./pDR('i)) (D 

where the free parameters (b, ai, a2, 03) are calculated using a 
Levenberg-Marquardt;^;-^ -minimization algorithm. A single optical 
depth value is applied to the three spectral components, but for most 
sources the solution is almost identical when two optical depth val- 
ues (one for the AGN and other for the PDR and HII components) 
are used, with an increased value due to the extra free parameter. 

t(A) is obtained from the Galactic Centre extinction law 
jChiar & T ielens 2006), while /agn is rep resented by the IRS spec- 
trum of the PG QSO 3C273 l lWu et al.l^009i) . /pdr by the ISO- 
CAM spectrum of a PDR in the reflection nebula NGC 7023 
( ICesarskvetai]|l996j) an d fmi by the ISOCAM spectrum of M17 
in the vicinity of OB stars jCesarskv et aljl996d) . Uncertainties in 
the fit parameters are estimated using Monte-Carlo simulations. 

We quantify the contribution, r, of each spectral component 
to the MIR spectrum as the ratio of its integrated luminosity to 
the total luminosity in the 5-15 i^im restframe range, tagn^ ''hii and 
'■pDR for the AGN, HII and PDR components. Figure |5] shows the 
distribution of the three spectral components in the sample. 

Only 235 of the 695 sources with available redshifts have IRS 
spectra covering the 5-15 |xm spectral range in its entirety. In many 
local galaxies, the IRS spectrum starts at 5.2-5.5 \im, while in 
higher redshift sources the long wavelength end of the spectrum 
is often shortwards of 15 [xm. Nevertheless, 379 spectra cover more 
than 95% of the 5-15 |,im range, and 595 show coverage ratios 
above 60%. 

To validate the results on sources with incomplete spectral 
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Figure 4. 5-15|-im spectra of the AGN, HII and PDR templates used in 
the specti'al decomposition. The silicate emission template used in the al- 
ternative spectral decomposition model for silicate emission sources is also 
shown. 

coverage, we have performed test runs of the decomposition algo- 
rithm on the subsample with full 5.5-15 |xm coverage, but restrict- 
ing the fitting range to 5-13.5 ^m, 5-12 |xm, and 5-10.5 |xm. For 
each spectral component, /, we define the change in the estimated 
share of the 5-15 |.im luminosity as: Ar, = r' - r,, where r'. is the 
value of r, obtained when fitting in the restricted range. 

The results (Table|4]l show that decreasing the long wavelength 
limit to 13.5 |,im has almost no effect, with typical Ar, of a few per- 
cent, and no significant bias. A further reduction of the long wave- 
length limit down to 12 significantly reduces the contribution of 
the HII template, by 7 percentage points on average, and increases 
the dispersion in all components by a factor ~ 3. Nevertheless, mov- 
ing the limit to 10.5 |xm does not increase the bias, which improves 
sligthly for the HII and PDR templates; but the dispersion further 
increases by another factor ~ 2, to almost 0.2 for HII, which ren- 
ders the decomposition results rather unreliable. Note that the HII 
component is by far the most affected by truncation of the spectrum 
at long wavelengths because of its steep slope. 

Arg y also shows high dispersion for a 12 |im limit, but most 
of the strong deviations originate in spectra with poor S/N which 
show large parameter uncertainties even with full spectral coverage 
(Figure[6j. There is no significant bias for Argj except for the test 
in the 5-10.5 \im spectral range. 

In order to obtain reliable measurements without sacrificing 
too many sources, we opt to require a lower limit of 60% coverage 
in the 5-10.5 i^im spectral range to perform spectral decomposition. 
In addition, we provide coverage ratios for all sources in the 
sample. 

Figure |7] shows 5sii versus the optical depth at 9.7 \im esti- 
mated from the b parameter (T97) for the 379 sources with >95% 
coverage in the 5-15|,im range. 

The tight correlation and the significantly larger uncertainties 
in 5sii compared to xgj are noteworthy. There appears to be a slope 
change at 5sii — 2, (T97 ~ 3) that might be caused by saturation 
of the silicate feature: since 5sji depends on the 9.7 |,im flux, it is 
difficult to obtain an accurate measurement when the flux level 
nears zero; also, emission in the 9.66 |xm S(3) line of molecular 
hydrogen (found in many ULIRGs) can increase the measured 9.7 
|,im flux. On the contrary, T9 7 takes advantage of the whole 5-15 
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Figure 5. Contribution of PDR and AGN components to the integrated 
5-15|.im restframe luminosity as derived from the spectral decomposition. 
Filled symbols indicate sources with full 5.5-15 |-im coverage in the IRS 
spectrum, while open symbols represent sources with only partial coverage. 
The amount of HII component can be calculated by the formula: rmi = 1 - 
('■agn+'Vdr)- 
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Figure 6. Comparison of decomposition rg 7 values obtained by trimming 
spectral coverage to 5-12|.im restframe (Y axis) versus rgj obtained using 
the full 5-15 |J.m range (X axis). The dotted line represents a linear fit to the 
data. 



|xm spectrum, and we can obtain a sensible estimate of the 9.7 
|xm optical depth from the wings of the silicate profile, even if the 
feature is saturated. 

When the silicate feature appears in emission, the best fit is 
obtained for a negative value of the opacity parameter b, since this 
produces a 10 \im bump with a shape alike that of the silicate emis- 
sion feature. This unphysical solution allows for a decent fit of sil- 
icate emission sources with no extra parameters in the model, and 



the associated negative T97 values correlate with 5sii values (Figure 
0. 

In order to check the validity of this approach, and to investi- 
gate the influence of a negative b value in the relative contribution 
of the three spectral components, we compare the results with an 
alternative spectral decomposition model for sources with the sili- 
cate feature in emission. In the alternative model, we fix the opacity 
parameter at = and include a fourth spectral template with the 
silicate emission profile of PG1351 (the source with the strongest 
silicate emission in the sample), obtained as the excess flux over 
the interpolated continuum between 8 and 13.5 ^im. Figure|4]shows 
the 5-15 [xm spectrum of the four components. The opacity needs 
to be fixed when using a silicate emission template to avoid degen- 
eracies, since an increase in the silicate emission component can be 
compensated by a simultaneous increase in the opacity parameter. 

To obtain a meassurement of the silicate strength we use the 
combined 9.7 ^m flux of the AGN, HII and PDR templates as the 
continuum estimate (F^^m), and the 9.7 |,im flux of the fourth (sili- 
cate) template as the peak flux estimate of the silicate feature (Fsu). 
The silicate strength is then calculated as S^,] = ln(l -t- Fsh/F^om)- 

Comparison of S^n with T9 7 from the original decomposition 
model shows a good correlation, with significantly less dispersion 
than the correlation between 5sii and T9 7 (Figure [7] inset plot). 
Comparison of ^Jii with 5sii also yields a good correlation, with 5s*i 
~ 0.91 Ssii. Changes in the relative contribution of spectral com- 
ponents, when adding the Silicate template, are small; in particular 
rpDR, a parameter we will use for diagnostics in ^6.21 is almost un- 
altered with a ArpDR averaging 0.0016 and a standard deviation of 
0.012. 

The alternative spectral decomposition model assumes the sil- 
icate feature peaks at ~9.8 i^tm, while there are instances in our 
sample where the peak of the emission is shifted towards longer 
wavelengths. These Ssn estimates should, therefore, be considered 
as an approximation to the real values. A more accurate meassure- 
ment would require the ability to shape and/or shift the silicate 
template in order to reproduce the diversity observed within the 
ATLAS sample. This and other issues related to the silicate feature 
in AGN will be addressed in a forthcoming work. 

The luminosity of the PDR component estimated in the fit 
(Lpdr) is a measurement of the total PAH luminosity in the galaxy. 
Comparison of Lp^R with the 6.2 |,im PAH luminosity (Figure[8j in- 
dicates the latter is responsible for roughly 10% of the t otal PAH lu- 
minos ity in the galaxy, in agreement with the results of lSmith et al.l 
( l2007h in a sample of local starbursts, and that this ratio does not 
significantly change in the higher luminosity sources. 

A comparison of EW6.2 with rpoR shows significant dispersion 
(Figure |9]l, mostly because of variations in the relative strength of 
PAH features. Unusually low 6.2 \im to total PAH luminosity ratios 
(L6.2/LPDR < 0.05) are found mainly in heavyly obscured ULIRGs, 
where an absorption band of water ice overlaps with the 6.2[im 
PAH emission, while very high ratios (L6 2/Lpdr > 0.2) appear in 
sources with peculiar MIR spectra that suggest unusual physical 
properties or artifacts in the data. 



6 MIR DIAGNOSTICS 
6.1 The fork diagram 

The "fork diagram", i.e. the strength of the 10 [^tm silicate feature, 
Ssii , versus the eq uivalent width of the 6.2 \im PAH feature, EWg 2, 
first discussed by ISpoon et alj ( l2007h . is now considered a stan- 
dard tool to distinguish between AGN- and starburst-dominated 
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Figure 7. Ssii versus Tg j for the sources in the sample with complete 5.5-15 
|.im spectral coverage. Linear fits for sources with Tgj < 3 and rgj > 3 are 
shown in dotted lines. The inset plot is a close up to the silicate emission 
region of the diagram where only sources with uncertainties in Ssu below 0.2 
are shown, for clarity. The open blue squares are the SsU estimates obtained 
from the alternative spectral decomposition model (see text), plotted as a 
function of Tg j . 
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Figure 9. EW6.2 versus rpDR for the sources in the sample with a S/N > 
2 detection of the 6.2 [im PAH feature. Symbols indicate the relative con- 
tribution of the 6.2 ^m to the integrated PAH luminosity as meassured by 
the L62/LpDR ratio: red squares indicate a ratio larger than 0.2, while green 
triangles mark ratios smaller than 0.05. 
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objects) or emission (type 1 objects) while the PAH features are 
less pronounced. 

A variant of the fork diagram can be obtained using the 11.3 
l^im PAH feature instead of the 6.2 |,im one (Figure [TOt. This ap- 
proach has the advantage of requiring less wavelength coverage and 
avoids the bias introduced in deeply obscured sources by water ice 
absorption in the wings of the 6.2 \im PAH feature (note that the 
11.3 i-im flux can be heavy ly depleted by silicate absorption, but 
since the silicate feature is very broad, the equivalent width is not 
significantly altered). In this variant of the diagram the starburst- 
dominated branch turns nearly vertical, but the topology remains 
unchanged. 

The parameters derived from the spectral decomposition allow 
us to build an alternative fork-like diagram, where tpdr substitutes 
the PAH EW and T9 7 takes the place of 5sii (Figure [Til. The advan- 
tage of this procedure is the results do not depend on the measure- 
ment of a single PAH band and the flux at the core of the silicate 
feature and the related caveats; instead, it uses the whole 5-15 |xm 
spectrum, making it less vulnerable to noise. 



Figure 8. Ratio between the 6.2 |.im PAH luminosity, as obtained by inte- 
gration in the 5.95-6.55 i-im band with corrections for flux lost in the wings 
of the profile (L52), and the integrated PAH luminosity as derived from the 
PDR component in specfi'al decomposition (Lpdr) as a function of for 
sources with a S/N > 2 detection of the 6.2 |,im PAH feature. The dotted line 
represents the best fit linear relation. 



sources (e.g. 'Farrah et al.l |2007|: IZakamska etal] |2008| ; IWu et al.l 
[2009; Willettet al. 2010). Starbursts tend to have large values of 
EW6.2, which decrease with increasing 5sii due to dep letion of the 
PAH band by absorption of adjacent water ice bands jSpoon et al.l 
l2002h and increased continuum emission, while in AGN-dominated 
galaxies silicates appear in shallower absorption (typically type 2 



6.2 MIR classification 

We use the spectral decomposition results to obtain a MIR classi- 
fication of the sources in the ATLAS sample. At a first stage, we 
separate AGN-dominated and starburst-dominated sources based 
on the fraction of PDR component, rpoR. 

We set the limit between AGN-dominated and starburst- 
dominated sources at rpoR = 0.15, which corresponds to EW62 ~ 
0.2 ^im according to the relation shown in Figure|9] 

Hernan-Caballero et al. (2009) showed that typical 
vL,,(5.5[.im)/LiRratios for AGN have values around 0.3 (their 
Figure 13), and starburst sources have mean L62/iiRratios of 0.01 
(their Table 9). Aproximating the 6.2 [xm continuum with the 5.5 
[xm luminosity would imply an equivalent width of EW6 2 ~ 0.2 
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Figure 10. Fork Diagram showing the 10 |.im silicate strength versus the 
11.3 \xra PAH equivalent width. Symbols indicate the classification of 
sources in the literature. Solid symbols indicate optical AGN classification 
as Seyfert 1 (blue circles), Seyfert 2 (yellow circles), intermediate Seyfert 
types (green circles), LINER (grey circles), quasar (blue boxes) and type 
2 quasar (yellow boxes); red open circles indicate Starburst or Hll-region 
type in the optical, while black crosses and triangles represent ULIRGs and 
submillimiter galaxies, respectively. Sources not classified, or with types 
that do not match any of the former ones are shown in black open squares. 











X » 








A 


® o 


X 

Xx'' 


□ 










» X 

> % " 

X X 






■ 




□ 


X 
X 

o* ® * X 
«x* *x 








■ 


X 
X 


X . x«|« 


■ 

o 


X 




1 

• • 

« 1 


• 


■ • «% 

X „ ^ 


X,* - 


• 
• 


■ 




' X •■ 
■ '0 


Jo° ■ 

□ 


s : 
■ • 
■ 


••i 


• 


^ ■ 













5.0 le-02 5.0 le-01 5.0 

("PDR 



Figure 11. Analog to the fork diagram obtained by representing the frac- 
tion of PDR component rpoR versus the optical depth at 9.7 |.im, rgj, as 
calculated with the spectral decomposition method. Symbols as in Fig. 1101 



|.im for a source with 50% starburst and 50% AGN contribution to 
the IR luminosity. 

Since many spectra are noisy and there is a significant number 
of sources near this threshold, we require the 1-cr confindence in- 
terval not to contain the threshold value in order to qualify as AGN- 
dominated or starburst-dominated. The remaining sources are clas- 
sified as "composites" if tpdr has S/N > 2 or "unknown" otherwise. 

In 93 sources the spectral coverage in the 5~15|,im range is too 
low to estimate rpuR, but the equivalent width of the 6.2 or 1 1.3 |.im 
PAH features is meassured. For these sources we use EWn.j = 0.2 
|.im (or EW52 = 0.2 |,im) as an alternative boundary between the 
starburst-dominated and AGN-dominated classes. 

With these criteria, 257 sources are classified as starburst- 
dominated in the MIR (MIR.SB), 348 as AGN-dominated 
(MIR_AGN) and 29 are too close to the threshold to decide, so they 
get classified as composites (MIR_COMP). The spectra of another 
49 sources are too noisy to tell, and in the remaining 12 the 6.2 and 
1 1.3 PAH features are both outside of the observed spectral range. 

The 348 AGN-dominated sources are further sub-classified 
according to their T97 meassurement into silicate emission 
(MIR.AGN1, 119 sources) or absorption (MIR.AGN2, 160 
sources), again requiring the 1-cr confidence interval in Tgj not to 
cross T9 7 = for a clean classification. Sources in which the l-cr 
interval contains Tgj = 0, but with uncertainties below 0. 1 are clas- 
sified as "flat MIR-spectrum AGN" (MIR_AGNx, 18 sources) to 
indicate that there is no strong emission or absorption at 10 |,im, 
and the remaining 5 1 sources retain the generic MIR_AGN classi- 
fication. 

Table [5] shows the distribution of sources in the sample ac- 
cording to their optical spectroscopic and MIR classifications. Val- 
ues indicate the number of sources with a given optical and MIR 
classification. Sources classified in the literature as both AGN and 
starburst in the optical are discarded for clarity, but the counts that 
would be obtained were they included, are shown in brackets. 

In most sources, the optical and MIR classifications are 
broadly consistent with each other (i.e., type 1 AGN usually show 
strong continuum and silicate emission, type 2 show also strong 
continuum and silicate absorption, and starbursts show strong PAH 
features on top of a fainter continuum), but there are exceptions. In 
particular, more than 1/4 of Seyfert Is and QSOs show significant 
silicate absorption, and many sources classified as AGN in the opti- 
cal are starburst dominated in the MIR. Apperture effects may be at 
least in part responsible for this, since the IRS spectrum of a galaxy 
nucleus include a far larger share of the emission from the galaxy 
compared to the corresponding optical spectrum. 

In the lower luminosity Seyferts, the host galaxy can con- 
tribute by a large fraction to the luminosity within the IRS slit, 
and composition of the (AGN-dominated) nuclear spectrum and 
(starburst-dominated) host galaxy spectrum can produce a silicate 
feature in absorption. This interpretation is supported by the strong 
PAH features observed in some Seyferts, that make them qualify as 
starburst-dominated in the MIR. Nevertheless, significant silicate 
absorption is also present in the spectrum of several QSOs, where 
the contribution from the host galaxy to the MIR luminosity is ex- 
pected to be small. Thus at least in some cases the silicate absorp- 
tion feature may be caused by obscuration of the AGN emission, 
be it in a dusty torus or in the host galaxy. 

On the other hand, there are also 15 Seyfert 2s and QS02s (1/6 
of the total) with the silicate in emission, as meassured by Tgj (12 
of them also confirmed by 5sii > 0). 

These apparent discrepancies, however, can not be regarded as 
evidence against the unified model of AGN; in fact, such behaviour 
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is predicted by both smooth and clumpy circumnuclear dust distri- 
bution models. Smooth model s predict silicate absorption in com- 
pact, high opacity face-on tori (Fritz et al. 2006), and also allow for 
silicate emission in type 2 AGN (e.g. Hatziminaoglou et al. 2008). 
Clumpy models require a high number of clouds along the equator 
and high optica l depth to produce silicate absorption in type 1 AGN 
jNenkova et al . 2008), while some parameter combi nations could 
give rise to silicate emission in type 2 AGN (e.g. iNikutta et al.l 
l2009l) . 

We therefore conclude that a silicate emission or absorption 
feature is not a reliable indicator, on its own, of the optical spectral 

type- 
There are also four LINERs with silicate emission, three of 
them show broad Balmer lines, and a strong silicate peak, while 
the remaining one is a radiogalaxy with a weaker silicate feature. 

Almost all optical starbursts are also classified as starbursts 
in the MIR, but the opposite is not true: roughly half of MIR star- 
bursts belong to some of the optical AGN classes. These sources are 
mostly nearby galaxies harbouring low luminosity AGN, that can 
be outshined by the galaxy in the large IRS aperture. When tak- 
ing into account only objects at z > 0. 1, from 98 MIR.SB sources 
54 are classified as starbursts in the optical, while 3 1 have LINER 
cores and only 13 are Seyferts or QSOs. 



starlight is not negligible compared to the AGN emission. Com- 
plementary templates are generated for all sources classified as 
ULIRG (h) in the literature (inespective of their optical spectral 
type) as well as the SMGs (;). The latter has a very limitted spec- 
tral coverage, as all the objects composing the SMG sample are 
high-redshift objects (with an average redshift (z) = 1.93). 

The MIR AGN average templates (j, k) also show interesting 
characteristics: the objects selected to have silicate in emission at 
10 [im have a flat spectrum with the silicate feature at 18 |.im clearly 
in emission, as well. The objects selected to have the silicate feature 
in absorption have a redder continuum and show no feature at 18 
[xm, be it in emission or in absorption. In a detailed comparative 
study between smooth and clumpy dust distributions, Fetre et al 
(in prep) show that such silicate features' behaviour is compatible 
with both clumpy and smooth dust distributions and hence we can 
not derive any constraints from the average templates. 

Interestingly, the average starburst template constructed from 
objects selected as such in the optical (g) shows more prominent 
PAH features and a redder continuum than the equivalent MIR tem- 
plate (/). This difference is most likely due to the different nature of 
the objects composing the two samples. The optical sample consists 
of low redshift starburst galaxies alone, while the MIR template has 
a higher average redshift, includes many ULIRGs as well as objects 
with significant AGN contribution. 



7 AVERAGE TEMPLATES 

A number of synthetic galaxy templates have been obtained from 
the average spectrum of populations matching certain criteria, 
shown in Table|6] 

We obtain average templates for Seyfert I, Seyfert I.X, 
Seyfert 2, LINER, QSO, QS02, starburst, ULIRGs and SMGs. 
ULIRGs with LINER or starburst cores are discarded when com- 
posing the LINER and starburst templates, that include only 
sources with yLv(7[xm) < lO'*'* erg s"'. Three additional templates 
are generated for sources classified in the MIR as dominated by 
starburst (MIR_SB), unobscured AGN (MIR_AGN1) and obscured 
AGN (MIR_AGN2). 

The individual spectra are normalized at 7 [im resframe and 
resampled to a common grid of wavelength values. The highest 
and lowest flux values in each wavelength bin are discarded, and 
the remaining ones are averaged with the arithmetic mean to obtain 
the template flux. The 12 templates obtained are shown in Figure 
II2I Shaded areas indicate the 1-cr dispersion of individual spectra, 
which is higher at longer wavelengths because of slope variations 
among sources of the same type. These templates are available as 
ASCII tables with columns for rest frame wavelength, normalized 
average flux density, l-cr dispersion in flux density and number of 
sources contributing to the average for each wavelength bin. 

The QSO template (e) is flat and almost featureless, while the 
type 2 QSOs (/) have on average a redder continuum and a some- 
what pronounced silicate feature in absorption. The equivalent with 
of the PAH features in the QSO and QS02 templates are much 
lower compared to those present in the average Syl (a) and Sy2 
(c) templates, in spite of th e high star formation rates deri ved from 
both the PAH featu res (e.g.lHeman-Caballero et"al]|2009t) and FIR 
luminosities (e.g. llfatziminaoglou et al.ll20IOl) . because the strong 
continuum emission from the AGN masks the underlying emission 
from the host galaxy. The LINER template {d) includes only nearby 
sources and its low PAHs equivalent widths are indicative of weak 
circumnuclear star formation. Its 5-9 |xm continuum is bluer than 
either the starburst or Seyfert templates because contribution from 



8 CONCLUSIONS 

We have collected a sample of 739 active and starburst galaxies 
with Spitzer/YRS spectra from the literature. The sample, though not 
statistical, comprises a variety of objects, namely type 1, 2 QSOs 
and Seyfert galaxies, intermediate Seyfert galaxies, ULIRGs, star- 
burst and sub-mm galaxies, with a large span of spectral properties 
such as PAH and silicate features, continua and distinct physical 
characteristics, and is the largest of its kind compiled to date. MIR 
properties were computed in a concise way for all objects and can 
thus be directly compared among the various types of objects. 

We find that the silicate strength, 5sii, a meassurement of the 
intensity of the 10 [im silicate feature that takes positive (negative) 
values when the feature is in emission (absorption), is well aproxi- 
mated by a "three band" color index, Cv-io-m, that is usually easier 
to obtain, except for sources with strong PAH emission that con- 
taminates the flux at 7 |,im, causing Cv-io-m to consistently lie short 
of ^sii. A more elaborate meassurement of dust obscuration, Tgj, 
obtained from the scaling factor applied to the extinction law in 
a spectral decomposition analysis, shows also a tight correlation 
with 5sii, but with different slopes for sources with saturated (T9 7 > 
3) and unsaturated (tj 7 < 3) silicate features. 

We find that the ratio of the PDR component to the total lumi- 
nosity of each object, integrated between 5 and 15 |,im restframe, 
rpDR, and T9 7 , both obtained from the spectral decomposition, pro- 
vide better AGN versus SB and obscuration diagnostics than the 
usual PAH EW and 5sii in sources with sufficient spectral coverage, 
as they do not depend on a single spectral feature and thus are less 
vulnerable to noise. Most QSOs and Seyfertl galaxies show T97 
values clustered around zero and tpdr values below 0. 15, while star- 
burst galaxies have large rpuR values and moderate to deep silicate 
absorption (higher in the starburst ULIRGs than in lower luminos- 
ity starbursts). Type 2 QSOs and Seyfert 2s exhibit low tpdr values 
and usually silicate absorption, sometimes strong, while LINERs 
can have almost any combination of parameters, due to their nature 
and definition. 
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Figure 12. Composite templates for subsamples of the ATLAS sample. In the first three rown, sources are selected based on their optical type and MIR 
luminosity, except for the ULIRG and SMG templates (see text). The bottom row shows average templates of the sources grouped according to the MIR 
classification. 



We then derive a simple classification of the sources, based on 
the values of rpoR and Tgj. Objects with rppR < 0. 15 are consid- 
ered to be AGN-dominated, while rpuR > 0. 15 indicates starburst- 
dominated sources, and the AGN-dominated subsample is further 
divided into silicate emitters and absorbers (tj 7 > and Tgj < 0, 
respectively). Comparison of this simple MIR classifications with 
the optical classifications in the literature is not straight forward, 
since the much wider apperture of the IRS slit with respect to op- 
tical spectroscopy allows for a larger fraction of the host galaxy in 
the MIR spectrum of resolved sources. Nevertheless, optical and 
MIR classifications are broadly consistent. Apperture effects may 
explain the larger fraction of starburst-dominated sources in the 



MIR relative to optical classifications, and the composition of nu- 
clear and host galaxy spectra might account for some of the type 
1 Seyfert galaxies that exhibit silicate absorption, even if in higher 
luminosity QSOs a source of obscuration (whether in a circumnu- 
clear dusty torus or in the host galaxy) is required to explain the 
cases of silicate absorption. Since up to 1/4 of the type 1 AGN 
show silicate absorption and 1/6 of type 2 AGN show silicate emis- 
sion, we conclude the silicate feature is not, on its own, a reliable 
indicator of the AGN optical type. 

All spectra in the sample, the average templates, as well 
as ancillary data, including (but not limited to) optical, IR 
and radio photometry, redshifts, and MIR spectral meassure- 
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ments are made available to the public at the ATLAS website: 
|http://www.denebola.org/atlas| 
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An atlas of mid-infrared spectra of galaxies 

Table 1. IRS observations log 



Sample IRS modules Spitzer Archive PID(s) Comments 



Brand et al. 2008, ApJ, 680, 1 19 


SL1,LL2.LL1 


15 


16 optically faint QSGs with X-ray detection 


Brandl et al. 2006 , ApJ, 653, 1 129 


SL.LL 


14 


22 nearby starbursts 


Buchanan et al. 2006, AJ, 132, 401 


SL2-LL1 (map) 


3269 


87 nearby 12um-selected Seyfert galaxies 


Cao et al. 2008, MNRAS, 390, 336 


SLJi 


various 


19 ulttaltiminous IR QSOs 


Dasyra et al. 2009, ApJ, 701, 1123 


SL1,LL2,H1* 


20629 


150 24[im selected IR-luminous galaxies 


Deo et al. 2007, ApJ, 671, 124 


SL4X 


3374 


12 Seyfert 1.8 and 1.9 galaxies 


Farrah et al. 2009, ApJ, 696, 2044 


SL1,LL2,LL1* 


30364 


16 optically faint 70 [xm sources 


Haas et al. 2005, A&A. 442, L39 


SL,LL 


3349 


7 quasars and 7 radio galaxies 


Hao et al. 2005. ApJ, 625, L75 


SL.LL 


14 


5 PG quasars 


Hernan-Caballero et al. 2009, MNRAS, 395, 1695 


SL,LL 


3640 


69 15um-selected IR-luminous galaxies 


Hiner et al. 2009 , ApJ, 706, 508 


SL,LL 


20083 


6 type 1 and 7 type 2 MIR selected quasars 


Imanishi et al. 2007, ApJ, 171, 72 


SL,LL 


105,2306,3187 


buried AGN in IR galaxies 


Imanishi et al. 2009, ApJ, 694, 751 


SLJi 


105,30407,20589 


20 buried AGN in IR galaxies 


Imanishi et al. 2010, ApJ, 709, 801 


SLJi 


50008,3187,30407 


17 nearby ULIRGs 


Uicy et al. 2007, ApJ, 669, 61 


SL1,LL2 


20083 


6 type 2 quasars 


LeipsM et al. 2009, ApJ, 701, 891 


SLi-L 


20525 


25 FR-I radio galaxies 


Lutz et al. 2008, ApJ, 684, 853 


LL 


30314 


12 z~2 mm-bright type 1 quasars 


Maiolino et al. 2007, A&A, 468, 979 


LLl 


20493 


4 high luminosity quasars 2<z<3.5 


Menindez-Delmestre et al. 2009, ApJ, 699, 667 


LL1,LL2* 


20081 


24 submm galaxies 


Murphy et al. 2009, ApJ, 698, 1380 


SL1,LL2,LL1 


20456 


22 galaxies in GOODS-N 


Netzer et al. 2007 ApJ, 666, 806 


SL,SH,LH 


3187 


23 Palomar Green QSOs 


PoUetta et al. 2008, ApJ, 675, 960 


SLl,LL2Jil 


various 


21 obscured AGNs 


Pope et al. 2008, ApJ, 675, 1 171 


Ii,2,LLl* 


20456,20081 


13 submm galaxies 


Shi et al. 2006, ApJ, 653, 127 


SLJX 


various 


68 AGN 


Siebenmorgen et al. 2005, A&A, 436, L5 


SL,H 


20231 


3C 249.1, 3C 351 


Sturm et al. 2005, ApJ, 629, L21 


SL,SH,LH 


3237 


NGC 3998 


Sturm et al. 2006, ApJ, 642, 8 1 


SL,LL* 


3223 


7 X-ray selected type 2 quasars 


Tommasin et al. 2008, ApJ, 676, 836 


SH,LH 


30291 


29 Seyfert galaxies 


Valiante et al. 2007, ApJ, 660, 1060 


SL.LL* 


3241 


9 submm galaxies 


Weedman et al. 2005, ApJ, 633, 706 


SL,LL,SH,LH 


14 


8 AGN 


Weedman et al. 2006a, ApJ, 651, 101 


SLl,LL2Jil* 


12,15 


24 optically faint 24 ftm selected AGN 


Weedman et al. 2006b, ApJ, 653, 101 


SLl,LL2Jil* 


15 


1 1 AGN and 9 starbursts 


Weedman et al. 2009, ApJ, 693, 370 


SLJX* 


40038,20128,20083,40539 


24 AGNs and some starbursts 


WiUett et al. 2010, ApJ, 713, 1393 


SL,LL,SH,LH 


30515,50591,105 


8 compact radio galaxies 


Wu et al. 2009, ApJ, 701,658 


SL,LL (map) 


various 


103 Seyfert galaxies 


Yan et al. 2007, ApJ. 658. 778 


SL1,LL2,LL1» 


3748 


52 24 pm seleted sources 


Zakamska et al. 2008. AJ, 136. 160 


SL1,LL2 


105,3163 


12QS02sfromSDSS 



*not all targets observed in all modules 



Table 2. monochromatic luminosity meassurement parameters 



band name 




bandwidth 






[|un] 


[|un] 




K 


2.2 


0.2 


29 


3.6 


3.6 


0.2 


102 


5.6 


5.55 


0.5 


485 


7 


7.0 


0.5 


630 


9 


9.0 


0.5 


656 


10 


10.0 


1.0 


624 


12 


12.0 


1.0 


565 


14 


14.0 


1.5 


508 


15 


15.0 


1.5 


489 


18 


18.0 


1.5 


430 


20 


20.0 


1.5 


397 


25 


25.0 


2.0 


352 



Table 3. PAH meassurement parameters 



6.2|rmPAH 11.3nmPAH 

central wavelength 6.25 11.3 

integration band 6.0-6.5 11.05-11.55 

continuum intervals 5.8-6.0,6.5-6.7 10.75-11.0,11.65-11.9 
assumed FWHM 0.2 0.2 
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Table 4. Bias in decomposition parameters caused by incomplete spectral coverage 





5.0-13.5|im 


5.0-12.0lun 


5.0-10.5|im 




mean 


a 


mean 




mean tr 


Al-AGN 


0.0029 


0.017 


0.0426 


0.051 


0.0590 0.118 


Arnil 


-0.0081 


0.034 


-0.0699 


0.098 


-0.0499 0.193 




0.0052 


0.029 


0.0272 


0.065 


-0.0092 0.096 


AT9.7 


-0.0005 


0.141 


0.0027 


0.560 


0.2723 0.820 



Table 5. Comparison of optical and MIR classifications 



MIR.AGNI MIRj^GNx MIR^GN2 MIRj^GN MIR.COMP M1R.SB 



Syl 


7 





3 








1(2) 


Syl.X 


26 


6 


18 


1 





20(24) 


Sy2 


7 


1 


25(26) 





1 


19(24) 


LINER 


4 





4 








9 


QSO 


61(63) 


8 


23(26) 


19 


2 


5(8) 


QS02 


8 


2 


31(34) 


7 


4 


5(7) 


Sbrst/HII 


0(2) 





5(12) 


1 


7 


74(89) 



Table 6. Composite spectral templates 



name 


Nsources 




<z> 






^max [l*m] 


comments 


Syl 


11 


0.002 


0.041 


0.205 


5.2 


24.6 


Seyfert 1 with vLv(7|im) < 10'^'^ erg s"^ 


Sylx 


72 


0.003 


0.091 


0.371 


5.0 


24.6 


intermediate Seyfert types (1.2, 1.5, 1.8, 1.9) 


Sy2 


53 


0.003 


0.045 


1.140 


5.2 


24.6 


Seyfert 2 with rLv(7nm) < lO*" erg s"' 


LINER 


16 


0.001 


0.034 


0.322 


5.2 


24.6 


LINER with yLv(7iim) < lO*" erg s"' 


QSO 


125 


0.020 


1.092 


3.355 


2.5 


24.6 


QSOl and Seyfert 1 with i'Lv(7nm) > lO"* erg s"' 


QS02 


65 


0.031 


1.062 


3.700 


3.6 


24.6 


QS02 and Seyfert 2 with vLv(7nm) > lO" erg s"' 


Sbrst 


16 


0.001 


0.091 


1.316 


5.2 


24.6 


Starburst or ffll with vLv(7nm) < lO" erg s"' 


ULIRG 


184 


0.018 


0.730 


2.704 


4.5 


24.6 


ULIRG (low and high redshift sources) 


SMG 


51 


0.557 


1.869 


3.350 


4.8 


12.0 


Submillimiter Galaxies 


MIR.AGN1 


119 


0.002 


0.455 


2.190 


4.0 


24.6 


MIR selected AGN with silicate emission 


M1R.AGN2 


160 


0.002 


0.549 


2.470 


4.5 


24.6 


MIR selected AGN with silicate absorption 


MIR.SB 


257 


0.001 


0.413 


2.000 


4.6 


24.6 


MIR selected starbursts 



